Vortices and phase transitions of interacting bosons in a rotating lattice ring 
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We investigate the ground state properties of the repulsive interacting bosons in an one dimen- 
sional rotating lattice ring, and reveal that the superfluid density of the system and the mass current 
it can carry in the rotating coordinate are periodic functions of the velocity of the lattice. In the 
weakly interacting limit, the vortices are generated at the critical velocities, and in the strongly 
interacting limit, the phase twist of the bosons induced by the rotating lattice manifests the super- 
fluid density and drive the system to undergo phase transitions between the superfluid and Mott 
insulator. 



The response of the superfluid to the external rotation 
is of great interest to the condensed matter physicists. 
Unlike the normal fluid, a superfluid can only follow the 
external rotation when the rotation frequency exceeds 
a critical value so that vortices can come into the sys- 
tem. The development of cold atom physics provides us 
with unique opportunity to investigate the critical ve- 
locity of superfluid through its highly controllability in 
experiments. For example, an optical lattice with de- 
sired lattice depth, geometry and rotation frequency can 
be applied to the cold atomic gas. Recent experiments 
at JILA have found the interesting pinning phenomenon 
and structural transition of the vortex lattice in an im- 
posed rotating two dimensional optical lattice^. A fun- 
damental question that arises naturally is, when and how 
do vortices come into the superfluid when a periodic ex- 
ternal potential is imposed? Moreover, it is known that, 
if the lattice barrier is increased, the superfluid will be- 
come weaker and weaker, and finally undergoes a phase 
transition from the superfluid to the Mott insulator at in- 
teger filling. This raises another question about the fate 
of the weak superfluid, can it survive under external ro- 
tation? Both of the two questions are not unique to cold 
atom physics, but also to the condensed matter physics, 
which involves the crystal environment. 

In this letter, we consider the ground state of the repul- 
sive interacting bosons in a rotating optical lattice ring. 
The ring geometry of the optical lattice has been pro- 
posed to be experimentally realizable by the interference 
of Laguerre-Gauss mode laser and a plane wave laserQ- 
Their proposal can be easily extended to produce a ro- 
tating lattice ring. If we open a detuning A between the 
frequencies of these two modes, the lattice ring will ro- 
tate at the speed of w oc , where R is the radius of the 
ring where the lattice sites locate and L is the number of 
sites. As shown in Fig.(P\), we assume the confinement 
of the system in z direction is quite weak, and in the di- 
lute limit, the behaviour of the atoms are sitting in the 
ground state in z direction . In the x — y plane, the lattice 
sites are located along a ring, and the confinement in the 
radial direction is strong enough to quench the motion of 
the atoms along this direction. So we will concentrate on 
the physics along the azimuthal direction, and the sys- 



tem can be simplified as a one dimensional moving lattice 
with periodic boundary condition. 

We would like to point out that our work is differ- 
ent from recent work in |2n4j i in which they studied the 
dynamical instability of a bosonic superfiuid with high 
crystal momentum, which is not the ground state of the 
system. Here we will consider the evolution of the ground 
state as a function of the velocity of the lattice. The 
main results of our work are summarized as following: 
(I) Both the superfluid density ps and mass current J 
in the rotating coordinate are periodic functions of the 
velocity v, i.e., J{v) — J{v + Ty) and Ps(v) — ps{v + T^), 
where Ty = where M is the mass of the bosons 

and a is the lattice spacing. (II) In the weakly interact- 
ing limit, vortices are generated at the critical velocities 
Vc = ''^MaiT^ through the Umklapp process, where I is an 
integer number. (Ill) In the strong interacting limit, be- 
fore the vortices come into the system, the rotating lattice 
drives the system from the superfiuid to Mott insulator. 
If the rotation is continuously increased, the superlfiuid 
phase comes back again with vortices 5]. 

Hamiltonian: To construct the proper basis for the 
many body hamiltonian, we start from the single particle 
hamiltonian in the moving coordinate, H = Hq — v ■ P, 
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FIG. 1: Left: The schematic figure for the system. Right 
top: the energy spectrum for L = 10. Right bottom: the 
interacting process scattering L particles initially in fc = 
state to k — ^ state. 
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where -ffo = 2M~'~^<»* ^® hamiltonian in the stationary 
coordinate, v is the velocity of the moving frame, and Viat 
is the lattice potential, Viat{r) = Viat{r + a). Hereafter, 
we use r to represent the coordinate in the ring. As 
H{r) = H{r+a), the solution of the Schrodinger equation 



Hil){r) — E->p{r) has the formula ipnkir) 



Unk{r) 



where Unkir) = Unk{r + a) and n, k is the band index 
and the crystal momentum respectively. To satisfy the 
periodic boundary condition in the ring, 'ipnk{r + La) = 
'4'nki'''), k can only take the discrete values, k — — 
0, ±1, ±2...± [L/2]. Substitute it back to the Schrodinger 
equation, we obtain 



(-iV - Mv + kf 
2M 



V — I Unk 



^nkUnk- (1) 



Compare it with the v — case, i.e., the equation for 
Unk in the stationary frame, we immediately know that 

2 

where the superscript denotes 



^nk 



^nk—Mv 



the physics observable in the stationary coordinate. In 
tight binding limit, e^'^j, = t„(l — cos/ca), we obtain 



e-nk = t„ (1 - cos (fe - Mv)a) 



(2) 



A typical energy spectrum within a single band as a 
function of v is shown in Fig.JQ). 

Now let us turn to the many body hamiltonian, 

H = j dr^\^^-v-P + Viat)^ + j dr^^m^^^, (3) 

where g — ^^^""^ ^ dzdijjUQ(z)hQ(uj) is the effective one 
dimensional interaction, a^c is the scattering length, 
uq{z) and ho{uj) are the ground state wave functions 
along z and radial direction. We can expand the sin- 
gle particle part of the Many body hamiltonian, K = 
J dr*t(^ - y . p + VJat)*, in the basis of Vnfc, ^' = 
J2nk '^nkdnk, whcrc a„fc is the annihilate operator for the 
single particle eigenstate i/'nfc- In this basis, K has a quite 
simple formula, K = J^nk^nkolik'^nk- In the coordinate 
space, it becomes 

K = - J2 inie''^biX,+c.c)+Y^it,,~^)blbm, 



n<i,j> 



where 6„i ^ ^J2k 



,-ikRi 



(4) 

a„fc is the annihilate opera- 



tor for the Wannier function Wni{r) = 



^-ikRi 



i'nkir) 



which is localized at site i and < i,j > denotes the near- 
est neighbor sites. We can also expand the interaction 
energy Eint = 9 J dr'if'' {r)'i>'' {r)'i> {r)'i> (r) in the basis of 
Wannier function. If we only keep the onsite interaction 
within a single band in the tight binding limit, we obtain 
the effective many body hamiltonian. 



where 9 = m = blbi, U = g J dr\Wi{r)\^ / . We 
have already dropped the constant term in the kinetic 
energy for simplicity. 

Small U limit: From Fig.|Q3), it is clear that, within 
each region ^^'^"^^^ < < ^2!d^llL^ single particle 
ground state i/'fc, ~ e*~'' carries a vortex with wilding 
number I, as the atom acquires a phase 2t:1 when it goes 
long the ring from r = to r = La . Turn to a many 
body system with N interacting particles, an interesting 
question arises immediately, how does the macroscopic 
wave function evolve from to -0 
the lattice velocity? 

To obtain the effective hamiltonian in the weakly in- 
teracting limit, we rewrite the hamiltonian in Eq.Q in 
the momentum space as. 



when we increase 



H 



E 



efc4«fc+y aiy,^^ak:,akJk^+k^-k3-k^-^, 
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(6) 

where s is an integer. Unlike the homogenous case, the 
difference between the total crystal momentum before 
and after the collision can be any reciprocal lattice vec- 
tor ^j^, which is well known as the Umklapp process 
in lattice. For simplicity, we concentrate on the region 
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In the weakly interacting limit, the 



ground state can be confined in the subspace composed 
by \m) = aJ'"al^'^~'"Vv/™!(^""^)!|0>, where al{a\) 
creates a particle in /cq = 0(fci = 2ti/L) state. Be- 
cause of the crystal momentum conservation, we notice 
{m'\H\m) ~ 5m,m'±L- To scatter L particle from fco state 
to fci state, we need incorporate the high order scatter- 
ing process, and the one of the lowest order is shown in 
Fig-lD^)- The effective hamiltonian is 

H = eoalao + eialai+goalaoa\ai+gi{a\^aQ +c.c), (7) 

where 50 = ^ and g^ = (^)^-Vn:=2 "'(^ - e,). 
One might assume the last term in Eq. is a high order 
process and should be quite weak, but as we have seen, 
the kinetic energy is degenerate at Vc- This coupling term 
thus plays a non-perturbational role. 

Rewrite the wave function as jG) — J2m Cm\N/2+mL), 
and project the Schrodinger equation to the Fock state, 
we obtain an effective Schrodinger equation for Cm, 



ECm = U{m)cn 



(Cm+1+C„ 



1), (8) 



where U^m) = ((eo — ei)mL + {g[ — go)ni'^L'^) and g[ — 
5i(f • At sufficiently large ffiling, i.e., {i^)^L > 

, it is straightforward to check that g[ will be greater 
than go. Consequently, the effective hamiltonian is an 
one dimensional harmonic oscillator and the ground state 
at Vc where cq = ci is a fragmental state 
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FIG. 2: Top: The superfluid fraction. Bottom: The mass 
current. The sohd, dotted and dash-dotted hne represent L — 



3, 6 and 10, where ^ = 0.9, and 
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keeps stationary in the laboratory frame. But when 
the lattice velocity approaches the first critical velocity 
-p^, more and more atoms occupy the state ■0fcij and 
both the amplitude of the mass current and superfluid 
fraction decreases. When the lattice velocity crosses Vc, 
most atoms will stay in the state ipki , and the many 
body wave function gradually changes to ■ The 
supefluid fraction comes back to 1 and a vortex enters 
the supefluid. If the velocity increases continuously, 
similar procedure repeats, vortices with different wilding 
number come into the system one by one. 

Large U limit: Wc will show that in this limit, the 
external rotation can drive the system from the superfluid 
phase to the Mott insulator. At large U limit, the number 
fluctuation in each lattice site is suppressed, and we can 
take the ansatz for the ground state. 



\G)^l[{a\n~l)+(]\n)+^\n+l}), 



(13) 



whose single particle density matrix has more than one 
macroscopic eigenvalues, 



(ajao) 
{a\ao) 



{a{ai) 











. 





(10) 



When the velocity of the lattice is away from the critical 
value Vc, Eq.(|SJ will be solved numerically to reveal how 
this fragmental state changes to the superfluid phases, 
which is reflected by evolution of the mass current and 
superfluid density. 

In the rotating frame, each atom in the eigenstate tpki 
carries the mass current Jk, — jf'^ki^ki m the rotating 
frame, and a straightforward calculation shows the aver- 
age mass current per particle is, 



J = 



Mta 



( (q£qo) 



sin(- 




(11) 



From the mass current, we can also define the superfluid 
fraction ris- As the normal part of bosons follow the 
lattice with a velocity u, the contribution to the non- 
zero mass current in the rotating frame totally comes 
from the superfluid fraction. When a vortex with wilding 
number I is present, the phase twist of the superfluid 
between two nearest neighbor sites is = ^g^, and the 
mass current can be written as J = nsMtasai{(j) — 0)/h, 
and consequently. 



hJ 



Mtasm{(j) - 9) ' 



(12) 



The numerical results are shown in Fig.lj^J. We see, at 
small lattice velocity, all the atoms stay at the state ipko 
and the superfluid fraction keeps as 1. The mass current 
in the rotating frame is negative, as the superfluid 



where jap + |/3p + = 1, and minimize the energy 
{G\H — ^J■N\G) for the grand canonical ensemble. The 
chemical potential /i is fixed by the total particle number. 
Because of the presence of the external rotation, the order 
parameter ipm = (bm) aquires a phase, and we can write 
it as ipm = i/Je'""^, where ip = |(&m)|. We have assumed a 
uniform phase twist in the ring. To satisfy the periodic 
boundary condition, (jj must satisfy {L — 1)0 — ~(f), i.e., 
(j> = = 0,±1,±2... ± [■§]. For simplicity, we wiU 

call the system is in the Ith branch if the phase (f> is 
We will see immediately that the phase 4> is not only 
the response of the superfluid to the external rotation, 
but also affects the amplitude of the superfluid density 
in turn. 

Apply the standard mean field procedure^, decou- 
pling the tunneling part of the hamiltonian as, 



we obtain the effective hamiltonian H 



(14) 



K 



T '>^m{nm " 1), whcrc 



K = J2 -2tcos(6'-(/))^(e*'^"6j„-fe-*"*'"6„0+2tcos( 



(15) 

Compared with the stationary case, it is clear that the 
tunneling is replaced by an effective one teS = t cos{9—4>). 
The factor cos(6' — 0) comes from the interference of two 
different ways which can generate an atom at site to, 
say, from site to — 1 and to -I- 1, and it also implies the 
periodicity of the all the properties of the system. On 
the other hand, it is well known that, in the stationary 
lattice, with the increase of U/t, the superfluid density 
will decrease. Especially, at integer filling, a phase tran- 
sition from superfluid to Mott insulator can occur. Con- 
sequently, the phase twist induced by the rotating lattice 
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FIG. 3: Top: The superfluid fraction. Bottom: The mass 
current. The sohd, dash-dotted and dashed line represent 
= 0.96, 0.98 and 0.99, where the number of sites is L = 10, 
and the particle number per site is n = 1000. 



reduces the effective tunneling and manifests the super- 
fluid density globally. This is quite different from the 
well known phenomenon, in a system without lattice the 
superfluid density can only be suppressed locally by rota- 
tion at the vortex core where the singularity of the phase 
occurs. 

If we start from a stationary lattice, and adiabatically 
increase the velocity of the lattice, the superfluid initially 
staying in the = state will vanish when 9 satisfy 
U = U°cos0, where U° = t{2n + 1 + ^(2n+ 1)^ - 1) 
is the critical interaction for the phase transition in the 
stationary coordinate. In another word, the transition 
happens when the velocity of the lattice reaches. 



h u 

— — arccos — . 
Ma t/P 



(16) 



From above equation, it is clear that vm can be much 
smaller than Vc = when U is sufficiently large, 

i.e., U > J7°cos^. In this large U limit, before the 
vortex emerge, the system has already entered the Mott 
insulator phase. If we continuously increase the velocity 
of the lattice, not surprisingly, the superfluid fluid phase 
will come back again, as the effective hamiltonian is a pe- 
riodic function of v. When -vm<v< jj^ -I- vm, 
the energy of Ith branch becomes the lowest one, and 
the superfluid state with a vortex of wilding number 
I is the ground state of the system. To quantitatively 
describe the system, we define the superfluid density as 
the usual way, ps = '■P^, and consequently, the superfluid 
fraction Us = Ps/n. According to Ea. (|12|l . the mass 



current can also be calculated. The results are shown 
in Fig.© by solving the hamiltonian in Ea. (|15l) . With 
the increase of the lattice velocity, / = 1, / ~ 2... 
branches come into the system one after another. Unlike 
the small U limit, the Mott insulator plays the role as 
an intermediate state between different superfluid phases. 

Conclusion and remarks: We have investigated 
the ground state properties of the repulsive interacting 
bosons in a rotating lattice ring. An effective hamiltonian 
in the rotating frame is deduced and various phenomena 
about the vortices generation and phase transitions of 
the bosons are found. As we mentioned at the beginning 
of this paper, the behaviour of the superfluid in a crystal 
environment is generally a fundamental problem not only 
for the cold atoms community, but also for the condensed 
matter physics. As an example, the recent experiments 
on supersohd in a cylinder [3,1a] have attracted a lot of at- 
tention. Although the microscopic physics there is quite 
different, as the crystal order there is built upon the in- 
teractions among the Helium atoms themselves, not from 
the imposed external potential, some interesting phenom- 
ena found there are surprisingly similar as what we found 
here, i.e., the unit circulation of lattice for the generation 
of the vortex in the system and the destruction of the su- 
perfluid density by the external rotation. The relation- 
ship between these two systems, as well as the underlying 
physics, remains to be investigated. 

We acknowledge helpful discussion with T.L Ho. This 
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